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We study the spin-charge separation in a Kondo-like model for an impurity with a spin and a charge (isospin) 
degree of freedom coupled to a single conduction channel (the "spin-charge" Kondo model). We show that 
the spin and charge Kondo effects can occur simultaneously at any coupling strength. In the continuum (wide- 
band or weak coupling) limit, the Kondo screening in each sector is independent, while at finite bandwidth 
and strong coupling the lattice effects lead to a renormalization of the effective Kondo exchange constants; 
nevertheless, universal spin and charge Kondo effects still occur. We find similar behavior in the two-impurity 
Anderson model with positive and negative electron-electron interaction and in the two-impurity Anderson- 
Holstein model with a single phonon mode. We comment on the applicability of such models to describe the 
conductance of deformable molecules with a local magnetic moment. 
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I. INTRODUCTION 

Recent measurements of the transport properties of indi- 
vidual molecules using scanning tunneling spectroscopy, me- 
chanically controllable break junctions, electromigration, and 
other methods have shown that the conductance strongly de- 
pends on electron-electron interaction as well as on the vibra- 
tional properties of electrode-molecule-electrode junctions, 
i.e. on the electron-phonon couplin g 1 ' 2,3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 . On the 
theoretical side, the role of phonons had been studied in 
the framework of the Holstei n 11 ' 12 and Anderson-Holstein 
model s - 13 ' 14 ' 15 ' 16 ' 17,18,19 ' 20 ' 21 ' 22 , where the charge couples lin- 
early to the displacement, and various models where the 
center-of-mass vibrations modulate the transparency of the 
tunneling barrier s - 23,24 ' 25 . 

For sufficiently strong electron-phonon coupling of the 
Holstein type and for suitably tuned electrostatic potential 
(gate voltage), the low-energy electron configurations may 
consist of the empty and doubly occupied molecular orbital, 
while the singly occupied states are only virtually excited 14 . 
When these conditions are fulfilled, a charge equivalent of the 
Kondo effect occur s 14 ' 15,16 ' 26 . The effective pseudospin de- 
gree of freedom is in this case the axial charge 27 or isospin 
I z = 1/2(Q — 1), where Q is the charge (level occupancy): 
empty orbital corresponds to isospin down, I z = —1/2, 
while the doubly occupied orbital corresponds to isospin up, 
I z = +1/2. The problem maps to the anisotropic charge 
Kondo model, where the impurity isospin couples to the 
isospin (charge and pairing) density of the conduction band 
via an isospin equivalent of the antiferromagnetic exchange 
interaction 15,26 ' 28 . At low temperature the isospin is screened 
by pairing fluctuations in the conduction band and the ground 
state is a complex many-particle Fermi liquid state which is 
an isospin singlet. A 7r/2 phase shift occurs for low-energy 
quasiparticle scattering and the molecule becomes fully con- 
ductive when the temperature is reduced below the corre- 
sponding charge Kondo temperature T§. In this respect, the 
charge Kondo effect is equivalent to the spin Kondo effect; 
the only difference are the interchanged roles of the isospin 
(charge and pairing) and spin degrees of freedom. We also re- 



mark that the charge Kondo screening of an isospin degree of 
freedom (electron pairing moment) is fundamentally different 
in its origin and its properties from the electrostatic screening 
of a point charge by the conduction electrons. 

The Kondo model is an effective one-dimensional quantum 
field theory since the magnetic impurity is assumed to cou- 
ple only to a one-dimensional continuum of conduction elec- 
tron states with s symmetry about the impurity sito^ 9 -. Low- 
dimensional field theories have unique properties due to topo- 
logical restrictions in reduced dimensionality; for example, 
fermions constrained to live on a ID line can scatter only for- 
wards and backwards. A notable effect in one-dimensional 
systems is the separation of electron spin and charge which 
had been intensively studied in Luttinger liquids 30 : fun- 
damental low-energy excitations are not charged spin- 1/2 
Fermi-liquid quasiparticles, but rather spin-1/2 neutral parti- 
cles (spinons) and charged spinless particles (holons). Such 
behavior has been found, for example, in one-dimensional 
solids such as SrCu02 31 and ballistic wires in GaAs/AlGaAs 
heterostructures 32 . 

Spin-charge separation also occurs in the Kondo prob- 
lem. Using bosonization techniques, conduction band fermion 
fields can be described in terms of spin-up and spin-down bo- 
son fields. These bosons correspond to the particle-hole ex- 
citations in the conduction band and they can be recombined 
to form separate spin and charge fields which are essentially 
independent, but subject to a gluing conditional which is 
the only remnant of the charge-1, spin-1/2 nature of physical 
fermion particles. In the single-impurity spin Kondo prob- 
lem, the impurity spin couples only to the spin field, while the 
charge field is decoupled 29,35 ' 36 . In this sense, the spin and 
charge degrees of freedom are separated. As commented in 
Ref.[33, the spin-charge separation allows the quantization of 
spin to persist even in the case of strong coupling of a quantum 
dot to the leads when the charge is no longer quantized. This 
explains why unitary conductance can be achieved in quantum 
dots at relatively high (Kondo) temperature, which is essen- 
tial for experimental observability of the Kondo effect in open 
quantum dots 38 . 

An interesting situation develops when one magnetic (spin) 



2 




Figure 1: (Color online) Schematic representations of models dis- 
cussed in the paper, a) Spin-charge Kondo model, b) Two-impurity 
Anderson model with positive and negative electron-electron interac- 
tion, c) Single-impurity Anderson-Holstein model, d) Two-impurity 
Anderson-Holstein model with electron repulsion in one orbital and 
electron-phonon coupling in the other. 



Kondo impurity and one charge (isospin) Kondo impurity cou- 
ple to the same conduction band, as schematically depicted in 
Fig- [It- Due to the spin-charge separation, two Kondo screen- 
ing cross-overs are expected to occur independently: the mag- 
netic moment on the first impurity will be screened at the spin 
Kondo temperature Tj| ~ exp(— 1/poJs) [po is the conduc- 
tion band density of states at the Fermi level, J5 is the spin 
Kondo antiferromagnetic exchange constant] and the pairing 
moment (isospin) on the second impurity will be screened at 
the charge Kondo temperature Tjf ~ exp(— 1/ poJc), where 
Jc is the charge Kondo exchange constant. 

The goal of this paper is to corroborate the prediction of 
separate and independent spin and charge Kondo effects. We 
study how this scenario breaks down in the narrow bandwidth 
limit (large Js/D and Jc/D) when lattice effects become 
important. We also study the related two-impurity Anderson 
model with one impurity with repulsive electron-electron in- 
teraction, U > 0, and one impurity with attractive interaction, 
U < 0, as we move away from the Kondo limit |f/|/7iT 3> 1 
by increasing the hybridization strength T (see Fig.HJ)). In the 
Anderson model, higher T implies stronger charge fluctuation 
on the magnetic (U > 0) site and stronger spin fluctuation on 
the isospin (U < 0) site. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the spin-charge Kondo model, describe its symme- 
tries and comment on the possible symmetry-breaking terms. 
In Sec. HIl] we present the numerical renormalization group 
(NRG) results for the thermodynamic quantities and study the 
degree of the spin-charge separation. In Sec. [IV] we study 
the model in the continuum (wide-band) limit using the non- 
Abelian bosonization to show how the spin and charge de- 
grees of freedom separate, while in Sec. [V] we study the op- 
posite limit of narrow bands. This naturally leads to the con- 



sideration of the strong-coupling fixed point in Sec. [VI] The 
spin-charge separation in two-impurity Anderson models is 
studied in Sec. lVIII while in Sec. lVIIll we comment on the ef- 
fects of the simultaneous spin and charge Kondo effect on the 
transport through molecular junctions. This is supplemented 
in Sec. [lX]by an explicit demonstration of the charge Kondo 
effect in the two-impurity Anderson-Holstein model. 

II. THE SPIN-CHARGE KONDO MODEL 

The conventional spin Kondo model (SKM) is22 

4^ + J s s(0) ■ S, (1) 

kfi 

where Js is the antiferromagnetic Kondo exchange constant, 
S ^ = \jp ^ 4/i ff w' c fcV = §E4°Wm' ( 2 ) 

c kk' [ifi' fifi' 

is the conduction-band spin density at the impurity site, 
where er are the Pauli matrices, cjL is creation operator for 
conduction-band electron with momentum k and spin p, N s 
is the number of states in the band, and /t = 1/N C ^ fc 
is the combination of the conduction band states which cou- 
ples to the impurity. S is the impurity spin- 1/2 operator. By 
analogy, the charge Kondo model (CKM) i&2& 

Hckm = ^2 £fcC L c fcM + J d{0) ■ I, (3) 
where Jc is the charge Kondo exchange constant, 

^°) = ^ E eL<w&'a' (4) 

c kk' olol' 

is the conduction-band isospin density at the impurity site and 
J is the impurity isospin-1/2 operator. Here ^ ka = {cL , Ck\ } 
is a Nambu spinoriSii! for the conduction band electrons, so 
that 

c kk' 

c kk' (5 ) 
= \(-ifUl+iflfl) 

**(°) = \jh (4 T c fe'T + Ai c -k'i -&kw) 

c kk' 

= ^(/ T t /T + / i t /l-l). 

It should be noted that the z-component of the isospin den- 
sity is the electron charge density measured with respect to 
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the half-filled band, while x and j/-components are the elec- 
tron pairing density. Contrary to the SKM, in CKM the impu- 
rity couples only to the charge sector of the conduction band, 
while the spin sector is decoupled (see Sec.lIVt. 

We now combine the two models and introduce the spin- 
charge Kondo model (SCKM) with one spin impurity and one 
isospin impurity, both located at the origin: 



#SCKM - ^2 efcC fcM Cfc ^ 

fe \l 



J s s(0) ■ S + J c i(0) • I. (6) 



and isospin operators 

k 



/ j c k-\ c k 



(7) 



An experimental realization of this model would consist of 
a single impurity with both charge and spin degrees of free- 
dom, such as a deformable molecule embedded between two 
electrodes. In one possible scenario, two different molecular 
orbitals are active. One is singly occupied and has a magnetic 
moment localized, for example, on a magnetic ion embedded 
in a molecule. The other orbital is as an extended molecular 
orbital which is strongly coupled to a local phonon mode. 

We remark that we have assumed the isospin part of 
the SCKM to be isotropic in spite of the fact that the 
Schrieffer- Wolff transformation applied to the Anderson- 
Holstein model in general yields an anisotropic effective ex- 
change interactio n 15,28 . The exchange anisotropy in spin- 1/2 
Kondo model is an irrelevant perturbation in the renormalisa- 
tion group sense 42 and we disregard it in the first part of the 
paper, since it plays no role in the context of the spin-charge 
separation. The anisotropy can be, however, important in the 
experimental realizations of the model since it enters the ex- 
pression for the charge Kondo temperature. We return to this 
point in Sec.lIXI 

The spin-charge Kondo model is distinct from the a — 
t Kondo model (also known as the compactified Kondo 
model ) 36 ' 43,44 ' 45 ' 46 , where a single spin degree of freedom is 
coupled via S ■ s(0) and S ■ i(0) terms to a single conduc- 
tion channel. In the compactified Kondo model, one takes ad- 
vantage of the spin-charge decoupling and uses the additional 
isospin degrees of freedom in one channel to mimic the spin 
degrees of freedom of the second channel in the two-channel 
Kondo model. Our model features two different degrees of 
freedom, has different symmetry and thus different properties. 
Nevertheless, our results on the effects of the finite band-width 
(lattice effects) on the spin-charge separation help understand 
why the compactified Kondo model is not fully equivalent to 
the two-channel Kondo model (see also Refs. l43ll4q) . 



A. Symmetries and symmetry-breaking terms 

We assume a linear dispersion relation = Dk for con- 
duction band electrons, where 2D is the band-width and the 
wave-number k ranges from —1 to 1. The band is half 
filled (the chemical potential is at [i = 0), so that we have 
particle-hole symmetry. The SCKM then has commuting 
SU(2) S pi n and SU(2)i sosp i n symmetries generated by the spin 



In the presence of the magnetic field h and of detuned elec- 
trostatic potential 5, we add the following perturbation terms 
to the Hamiltonian: 



H' = hS z + 5(Q-l) = hS z + S (21*). 



(8) 



Magnetic field breaks the SU(2) sp i n symmetry, while detuned 
electric potential breaks the SU(2)i SOSp i n symmetry. Both 
perturbations are relevant: strong magnetic field h > Tj| 
quenches the spin Kondo effect, while strong potential 6 > 
T2 quenches the charge Kondo effect. 



III. NUMERICAL RESULTS 

We studied the model using the well-established numerical 
renormalization group (NRG) technique 47,48 . We used an im- 
plementation of the NRG where the conservation of spin and 
isospin is explicitly taken into account. 

We should emphasize the assumption of linear dispersion 
for conduction band electrons. It is known that in one- 
dimensional models spin and charge degrees of freedom truly 
separate only if the dispersion is exactly linear. The separa- 
tion actually extends to the region of non-linear dispersion, 
however the spin and charge excitations acquire a finite life- 
tim o 30 ' 49 ' 50 ' 51 . By performing the calculations for a range of 
the discretization parameters A 48,52 , we have verified that the 
results converge rapidly to the continuous band limit (A — > 1) 
and that the logarithmic discretization of the conduction band 
does not spoil the linearity. In this work we thus neglect the 
effects of non-linear conduction band dispersion, but we note 
that they could in fact also be studied using suitably adapted 
NRGcode 5 ^. 

We computed the impurity contributions to the entropy and 
to the spin and charge susceptibilities. These thermodynamic 
quantities are defined as^ 4 . 



Xspin(7 n ) 



(E -F) (E- F) 
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k B T 
1 



((Si) - (S*)o) , 



(9) 



Xchargc(T) = — «/ 2 ) - (7|) ) , 



where the subscript refers to the situation when no impu- 
rities are present, E = (H) = Tr (Her H / {k ' BT ^), F = 
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Figure 2: (Color online) Impurity susceptibilities and entropy for 
the spin-charge Kondo model. Filled symbols are a fit to the ex- 
act Bethe-Ansatz results for susceptibility. From this fit we de- 
duce Tk = Tk — 2.7 x 10~ 7 D, in approximate agreement with 
T K = £>^oJexp(-l/pJ) = 4.4 x 1(T 7 £>. 
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Figure 3: (Color online) a) Spin susceptibilities in the SCKM model 
with Jc = Js (solid lines) compared with spin susceptibilities in 
the corresponding SKM with the same Kondo exchange constant Js 
(dashed lines), b) Effective exchange constant JseB as a function of 
Js for the case Jc = Js (red, circles), compared with the single- 
impurity results (black, squares). Blue straight line is a fit to the 
single-impurity results and serves as a guide to the eye. c) Effective 
exchange constant JseS as a function of the charge Kondo exchange 
constant Jc while keeping the bare spin exchange constant Js con- 
stant at pJs = 0.1 (black, circles) or pJs = 0.05 (red, squares). 



— ksT In Tr (e _H / fefsT ) , g is the electronic gyromagnetic fac- 
tor, hb the Bohr magneton and kg the Boltzmann's constant. 

We calculated these quantities as a function of the tempera- 
ture for a number of choices of Js and Jc, see Fig. [2] We used 
discretization parameters A = 2,3,4 and six different val- 
ues of parameter z with the discretization scheme described in 
Ref.|53; the results nearly overlap for all three A. At each iter- 
ation we kept states up to an energy cut-off of at least 15Tjv, 
where Tjy is the energy scale at the A^-th iteration. Even-odd 
effects are removed by averaging over two consecutive itera- 
tions. 

As expected, we observe Kondo screening of both spin 
and charge degrees of freedom. Kondo effects in each sec- 
tor appear fully independent in the sense that the susceptibil- 
ity curves Xspin and Xchargo follow the universal S = 1/2 
Kondo forms and each is characterized by a single parame- 
ter, the Kondo temperature or Tj^, respectively. Even for 
Js = Jc, when the screening in each sector occurs at the 
same temperature, there is no competition between the two 
sectors and the curves agree perfectly with the exact Bethe- 
Ansatz results for the conventional S = 1/2 Kondo model 
(filled squares in Fig. [2]). 

On a lattice, spin and isospin cannot simultaneously be fi- 
nite at one sit o 36 ' 46 . We therefore expect that for large Jc and 
Js, when the Kondo temperatures are high and the Kondo 



screening clouds are small, there should be some kind of mu- 
tual disturbance. We performed a series of calculations for 
increasing parameter J = Js = Jc an d compared the spin 
susceptibility with that of a single-impurity spin Kondo model 
with the same Js, see Fig. [3^. We find that the universality 
is maintained: for any J, the susceptibility still follows the 
universal Kondo curve. The Kondo temperature is, however, 
reduced in the spin-charge Kondo model and it is no longer 
accurately described by the familiar expression 

T K ~ Dy/pJe'^. (10) 

A convenient quantity to study this renormalization effect are 
effective parameters Js e f£ an d JceS, which can be obtained 
from respective Kondo temperatures and Tjf using the in- 
verse of Eq. ( [Tol l: 

2 

pJcS ~ w[2/(T K /Dfy 

where W(z) is the product logarithm function (the solution 
x of the equation z — x exp x). In Fig. \3]p we plot the Js e fi 
which corresponds to the susceptibility curves from Fig. [3^. 
For Jcl) (i.e. in the wide-band limit), there is no renormal- 
ization and J c ff k J. For J < D we observe a systematic re- 
duction of the effective exchange constant for increasing bare 
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Js- This can be understood in terms of the "stiffness" of the 
Kondo clouds in spin and charge sector. With increasing J, 
the screening clouds are made to occupy smaller and smaller 
spatial extent around the impurity. Surprisingly, this does not 
lead to a collapse of the Kondo effects, as might be expected. 
Instead, the Kondo clouds spread out more than they do in the 
single impurity case to compensate for the conduction elec- 
tron spin and isospin density loss due to the lattice effects. 

In Fig. [3}; we show the variation of Js c s for a constant bare 
Js as the charge Kondo exchange Jc is increased from 0. The 
renormalization of the spin exchange constant is weaker for 
smaller bare Js, when the spin Kondo cloud is more spread 
out, which confirms our interpretation in terms of interfering 
Kondo clouds. 



IV. CONTINUUM (WIDE-BAND) LIMIT AND 
BOSONIZATION 



The symmetry of the problem and the separation of the spin 
and charge sectors is directly exhibited in the non-Abelian 
bosonization approach of the boundary conformal field theory 
(CFT)^. The 

essence of this approach is to represent the 
electrons as independent bosonic fields carrying the isospin 
SU(2) and spin SU(2) degrees of freedom 45 i 56 i 57 . Introduc- 
ing the left-moving field operators ip^fx), we write the spin- 
charge Kondo Hamiltonian in real space as 

iv Fsr f°° t , gVy(g) , 

Hsckm = — ^ y : ^(x)-^- : dx 
+ v F p(x s J s (0)-S + X c J C (0)-l) , 



(12) 



where vp is the Fermi velocity, p is the density of states and 
J s (x) and J c (x) are spin and charge (isospin) "currents' 
defined as 



J S ( X ) = ■ ^li^o-^uix) :, 

jC ( X ) = : ZiWoVcxP&ix) ■ 



(13) 



where ^(x) = ^ij^(x),ipi(x)^ is the real-space Nambu 

spinor. Normal ordering (double dots) has been introduced 
to remove divergences due to filled electron levels below the 
Fermi level. Note also that in this section the spin and charge 
Kondo coupling constants are Xs and Ac. 

The Sugawara form 58 equivalent to the Hamiltonian 
Eq. (Ell is 5 ££ 



Hi 



SCKM 



n n 

^(xsYjn-S + XcYjn- 1 ) 

\ n n / 



C % jC 



Here and are the Fourier modes of the real-space cur- 
rents J (x) and J (x) in a finite system of length 21. They 
each satisfy SU(2)i Kac-Moody commutation relations 

[J^ C ' a , Jm C,b ] — ^abcJn+m + ^ab^n+m,0 (15) 

where e abc is the antisymmetric tensor 56 . The modes from 
different sectors commute: 



jS,a jO,bl _ q 



(16) 



This commutation relation embodies the (trivial) spin-charge 
separation of the free electrons (Xs = Xc = 0). At a special 
value As = 1 /3, we can introduce a new current J s = J s + 
S which satisfies the same Kac-Moody commutation relations 
as the old currents. The spin part of the Hamiltonian then 
becomes (up to a constant term) 



(S) 

SCKM 



nv F x ^ 1 
3 



•J —71. 



J, 



S . 



(17) 



from which the spin impurity S has disappeared (it was "ab- 
sorbed" by the conduction band). The charge sector remains 
unaffected by this change. By analogy, at a special value 
Ac = 1/3 we introduce J c — J c + I and write the charge 
part of the Hamiltonian as 



H 



(C) 

SCKM 



TTVp 



El 



•~> — n >Jn 



(18) 



The special values As = 1/3, Ac = 1/3 are identified with 
the strong coupling fixed point of the problem 29 . Note that 



J, 



'S,a qC,b 



0. 



(19) 



(14) 



therefore the addition of the impurities does not break the 
spin-charge separation. 

Even though the spin Kondo effect occurs in the spin sector, 
without involving the charge sector, and the charge Kondo ef- 
fect occurs in the charge sector, without involving the spin 
sector, the spin and charge degrees of freedom are not en- 
tirely decoupled; they are constrained by the gluing condition. 
The gluing condition declares which combinations of quan- 
tum numbers are allowed taking into account the charge- 1 
spin- 1/2 nature of physical particles - electrons. In the present 
context the gluing condition depends on the boundary condi- 
tions (b. c.) imposed on the field ip(x). 

We first consider the case of anti-periodic b. c, t/j(l) = 
—ip(—l). We can obtain half-integer spin only with an odd 
number of electrons (i.e. for half-integer isospin) 29 . Therefore 
2I Z and 2S Z must have the same parity; this is the gluing 
condition. There are thus two Kac-Moody conformal towers 
with highest-weight states having (I,S) = (1/2,1/2) and 
(I, S) = (0, 0), respectively. 

For periodic b.c. ip(l) = tp(—l), and keeping in mind 
that the z-component of the isospin is defined with respect 
to half filling, we obtain half-integer spin for integer isospin 
and integer spin for half-integer isospin; 2I Z and 2S Z must 
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then have different parity. There are two conformal towers, 
(I,S) = (1/2,0) and (I, S) = (0,1/2). Note that changing 
the b. c. from periodic to anti-periodic (or vice versa) amounts 
to imposing a phase shift of tt/2 on the wave function 55 . 

In the single impurity Kondo model, the finite-size spec- 
trum of the strong coupling fixed point is obtained by a fu- 
sion in the spin secto r 29 ! 55 . This means that the isospin sec- 
tor remains intact, while the spin quantum number changes 
as S -> 1/2 - S. As a consequence (1/2,1/2) -> (1/2,0) 
and (0, 0) — > (0, 1/2), i.e. the gluing conditions change from 
those for the anti-periodic b. c. to those for periodic b. c, and 
vice versa, which is equivalent to the tt/2 phase shift. 

We now generalize this fusion rule to the case of both (spin 
and isospin) impurities. The absorption of the spin impurity 
will not affect the isospin sector (/ — > /), however S — > 1/2 — 
S. Similarly, the absorption of the isospin impurity does not 
affect the spin sector S — > S, while I —* 1/2 — I. The fusion 
rule is thus S -> 1/2 - S and / -»■ 1/2-1. Therefore 
(0, 0) -f (1/2, 1 /2) and (1/2, 1 /2) -v (0, 0) for periodic b.c. 
and (1/2,0) -> (0,1/2) and (0,1/2) -» (1/2,0) for anti- 
periodic b.c: the boundary conditions remain the same, only 
the conformal towers are permuted. There is no phase shift 
in this case. Alternatively, we can consider the fusion in each 
sector to give a tt/2 phase shift. This gives a total phase shift 
of 7r which, however, is equivalent to zero phase shift since 
phase shifts are defined modulo tt. The finite-size spectrum is 
that of free fermions, therefore the fixed point corresponds to 
a Fermi liquid. 



V. ZERO-BANDWIDTH AND NARROW-BAND LIMITS 

For very large Jc and J5, SCKM reduces to the zero- 
bandwidth limit where the band is effectively described as a 
single orbital. In this three-site problem, the ground state for 
Js > Jc is a spin singlet + a free isospin, while for Jc > Js 
the ground state is an isospin singlet + a free spin. Spin and 
isospin are in direct competition since a single band "orbital" 
can either behave as a spin degree of freedom or as an isospin 
degree of freedom, but not both at the same time. Only one of 
the two possible Kondo ground states can be realized. 

At finite bandwidth additional degrees of freedom become 
available and the spin and charge Kondo effects can occur si- 
multaneously. In the next approximation, we therefore take 
two lattice sites to mimic the conduction band. This "narrow- 
band limit" is the minimal model for a band with indepen- 
dent spin and isospin degrees of freedom. We use -ffband = 
— i^ucj^cfyj + H.c, where t <C Jc,Js- In this approxi- 
mation the ground state is spin-singlet isospin-singlet for any 
values of Jc and Js, as expected. The ground state energy is 
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Figure 4: (Color online) Schematic representation of the strong cou- 
pling fixed point Hamiltonian of the spin-charge Kondo model for a) 
Js > Jc and b) Jc > Js- F° r Jc = Js, the fixed point corre- 
sponds to the symmetric combination of a) and b). 



which, for large J s /t, Jc /t, equals 



E, 



GS 



-f max{ J c , J s } ~ I \j c t -j s \ if l J c - Js\^>t, 
y -\J-t if\J c -Js\<£t. 

(21) 

In the second line J ~ Jc ~ Js- The fixed point Hamiltonian 
consists of the impurity sites strongly coupled to the first two 
sites of the Wilson chain, which thereby become decoupled 
from the remainder of the chain. 



(20) 



VI. FIXED POINTS 



In the single-impurity Kondo model with one spin (isospin) 
impurity, the strong coupling fixed point corresponds to the 
impurity site and the first site of the Wilson chain being tightly 
bound into a spin (isospin) singlet, which results in the decou- 
pling of the first site of the Wilson chain from the res t 47,48 . As 
a result, the low-energy electrons experience a 7r/2 scattering 
phase shift. As shown in the previous section, in the strong 
coupling fixed point of the SCKM two sites are removed from 
the chain (see Fig.|4|. This corresponds to a zero phase shift 
for quasiparticles, in agreement with the CFT analysis. 

The SU(2) S pi n x SU(2); sosp i n ~ 0(4) symmetry with zero 
phase shift of the SCKM model should be contrasted with the 
0(3) x 0(1) symmetry of the a — r model at its NFL fixed 
point or the 0(4) symmetry with 7r/2 phase shift at its FL 
fixed point when an anisotropy between spin and isospin sec- 



tors exist: 



,43.45 
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VII. SPIN-CHARGE SEPARATION IN THE 
TWO-IMPURITY ANDERSON MODEL 

The Kondo model is an effective low-temperature theory of 
a magnetic impurity. If the physical reality is modelled more 
accurately, the impurity is described using the single impu- 
rity Anderson model (SIAM) 59 and the Kondo model arises 
only as a low-temperature effective theory of SIAM after 
the charge-fluctuations are projected out using the Schrieffer- 
Wolff transformation 60 . In SIAM, charge fluctuations persist 
at low temperatures and, in fact, they must be present since 
they provide the mechanism by which the impurity spin can 
flip. Nevertheless, in the strong Kondo regime, U/irT 3> 1 
(where U is the electron repulsion and T the hybridization 
strength), the Anderson impurity still couples predominantly 
to the spin sector of the conduction band and only weakly to 
the charge sector. In fact, this problem has to be considered 
from the renormalization group point of view. At high energy 
(temperature), when the system is in the free-orbital regime, 
the impurity indeed couples to both spin and charge sectors. 
At low energy, when the system is in the local moment regime, 
the coupling to the charge sector is frozen out. This implies 
that the spin and charge separate only for low energy elec- 
trons, while the two degrees of freedom are "entangled" for 
high energy electrons. Due to the energy-scale separation in 
quantum impurity models 47 , the lack of the spin-charge sepa- 
ration at high energies does not preclude the spin-charge sepa- 
ration at low energies (i.e. at low temperatures). Alternatively, 
this can be phrased in terms of the separation of time scales 54 : 
the duration of fluctuations (spin flips) is tjj ~ h/U while the 
"magnetic" time scale (roughly equivalent to the mean time 
between successive spin flips) is tk ~ ^/Tk- In the Kondo 
limit, Tk <C U and therefore tk 3> W. 

In the Anderson model with negative U the spin and 
charge sectors are interchanged 2 ^. For temperatures below 
\U\, a pairing moment develops in place of the magnetic mo- 
ment. The impurity state can be flipped from zero-occupancy 
(isospin down) to double-occupancy (isospin up) by coupling 
to the conduction band. The single-occupancy states can be 
projected out using a suitably generalized Schrieffer- Wolff 
transformation and the effective model is the charge (isospin) 
Kondo model^ 6 .. 

Generalizing the spin-charge Kondo model, we now study 
the two-impurity Anderson model (2IAM) described by the 
Hamiltonian 

H = ekC^Ckc + ^(ni - l) 2 + 5i(n : - 1) 

ka 

+ y(™ 2 -l) 2 + <5 2 (n 2 -l) 

+ J2 V£ (dLcfa, + H - c -) + E V k { d tc k « + H.c.) . 

ka ka 

(22) 

Here rii = ^ CT <$ iG di a is the electron occupancy of impurity 
i, U\ > is the electron-electron (e-e) repulsion on impurity 
1, U2 < is the e-e attraction on impurity 2. We assume a 
constant hybridization strength T = 7rpo|Vfe F | 2 ; this permits 
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Figure 5: (Color online) a) Spin susceptibility for 21 AM (full 
lines). We choose U\ = -U 2 = U and Ti = T 2 = T. 
For comparison we plot results of the SIAM with the same T and 
U (dashed lines). Parameters are U/D = 0.01 and F/U = 
0.02,0.03,0.04,0.05,0.1,0.2 and 0.5 (left to right), b) Charge 
fluctuations and correlations in 2IAM for increasing hybridization 
strength T. 



comparison with the Kondo model which, in a similar man- 
ner, has a constant density of state po and a constant exchange 
constant Jk (this corresponds to approximation Vk = Vk F ). 
Parameters Si = e, + C/j/2 measure the departure from the 
particle-hole (p-h) symmetric point at 6% = #2 = 0. The 
model is schematically depicted in Fig. QJ. In the p-h sym- 
metric point, 2IAM has SU(2) sp j n x SU(2)j sosp i n symmetry 
likeSCKM. 

In Fig. [5] we show the spin susceptibility and expectation 
values of the local charge fluctuations Anf = (n^ — (n>i)) 2 , 
inter-impurity charge fluctuations Ani?T,2 = (ni — {nx))(ri2 — 
(112)), and the inter-impurity spin correlations Si ■ S2 as a 
function of the hybridization strength V. 

For small T, i.e. for small equivalent Kondo exchange cou- 
plings 

8Ti T 8r 2 



Js 



Jc 



n\U 2 y 



(23) 



the first impurity behaves at low temperatures as the SIAM 
with the same Ui,Ti, see Fig. [5^. The good agreement at 
low temperatures demonstrates that spin and charge degrees 
of freedom in the 2IAM effectively separate in spite of charge 
fluctuation on the impurity 1 and spin fluctuations on the im- 
purity 2, as discussed above for the case of a single Anderson 
impurity. At temperatures above the moment formation tem- 
perature T* k 1/5C/ 48 , the magnetic susceptibility is larger 
than that of the SIAM model by a factor of 2 as both impuri- 
ties are then in the free orbital regime with equally probable 
configurations |0), | T), I I), I ?!>• 

For increasing T/U, the renormalization effect due to the 
presence of impurity 2 becomes noticeable and the Kondo 
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Figure 6: (Color online) Numerical renormalization group eigen- 
value flows for the case when a) isospin impurity only, b) spin im- 
purity only, c) both impurities are present. Ui = U = 0.01, 
r, 1, 2/U — 0.05 (unless turned to zero in a and b), U2 = —3(7. 
The levels are labeled by the total isospin and spin quantum num- 
bers, (7,5"). 




Figure 7: (Color online) Impurity susceptibilities and entropy for the 
two-impurity Anderson model. Impurity 1 has Ui — U = 0.01, 
T/U = 0.05. Impurity 2 has U 2 = -37/ = -0.03 and equal 
hybridization T. T^/D ~ 1.1 10" 6 , T%/D ~ 2.9 10" 13 . Note 
that the magnetic susceptibility curves nearly overlap for all values 
of 5. 



temperature is lower than in the corresponding SIAM model. 
This can again be explained in terms of Kondo cloud "stiff- 
ness", see Sec. [TTTJ For very large T/U we observe an inter- 
esting saturation effect. While for U/irT < 1 the Kondo effect 
in SIAM collapses, charge fluctuations rise to the maximum 
value of An 2 = 0.5 and the susceptibility drops to zero at 
T ~ r 48 , in 21 AM the charge fluctuations on impurity 1 are 
constrained due to the presence of impurity 2 (see the slow rise 
of the Anj curve in Fig. |5j>) and, similarly, spin fluctuations 
on impurity 2 are constrained due to the presence of impu- 
rity 1 (see An| curve). The temperature dependence of the 
susceptibility is in this case approximately T\{T) ~ logT 
in a temperature interval of four orders of magnitude, with 
the low-temperature tail still having the form of the universal 
Kondo curve. 

In Fig. [6] we show NRG eigenvalue flows for: a) single 
77 < impurity, b) single 77 > impurity and c) the 2IAM 
with both impurities. It clearly shows the difference between 
SIAM and 2IAM: for negative- 77 SIAM the ground state is 
an isospin singlet and the system is a Fermi liquid with a 
<5 qp = 7r/2 phase shift; for positive- 7/ SIAM the ground 
state is a spin singlet and the system is a Fermi liquid with 
a <5q P = 7r/2 phase shift; finally, in the 2IAM the ground state 
is a spin-singlet isospin-singlet (for odd N; for even TV that 



we show, the finite-size ground state consists of degenerate 
/ = 1/2, 5" = 0, I — 0, S = 1/2 states) and the system is a 
Fermi liquid with zero phase shift as in the SCKM model, see 
Sec. ED 

The departure from the p-h symmetric point is a marginal 
perturbation for the positive-77 impurity and, as long as \8± | < 
77/2, the spin Kondo effect persists 48 . The p-h symmetry 
breaking is, however, a relevant perturbation for the negative- 
77 impurity, where it plays the same role as a magnetic field 
for the positive-7_7 case (see Sec. Ill Al l. If \S 2 \ > T%, the 
charge Kondo effect is washed out. This is illustrated in Fig. [7] 
where we plot the thermodynamic quantities for a range of 
5 = 61 = 62- It should be noticed that the spin susceptibil- 
ity is hardly affected; the curves nearly overlap. The effect of 
a magnetic field is analogous: it washes out the spin Kondo 
effect, but it is marginal for the charge Kondo effect. 

In Fig.|7]we observe some high temperature features in the 
thermodynamic properties at T ~ U\, 77 2 , e.g., a dip in the 
spin susceptibility and a peak in the charge susceptibility. This 
is a common feature of the 2IAM models with U\ ^ —U-2- 
In a single impurity model with U\ > 0, the spin suscep- 
tibility increases at the local moment formation temperature 
T* = 1/577! from its free orbital value of 1/8 to the local 
moment value of 1 /4, while the charge susceptibility sharply 
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drops to zero as the charge fluctuations are frozen out. In a 
single impurity model with U2 < the same scenario oc- 
curs at T 2 * — 1/5 1 £^2 1 with the spin and isospin sectors ex- 
changed. Therefore, in 2IAM with U\ < | C/2 1 , as in this 
case, we observe a dip in the spin susceptibility and a peak 
in the charge susceptibility, since the pairing moment begins 
to form at a higher temperature than the magnetic moment. 
For U\ > I U2 1 , there is a peak in the spin susceptibility and a 
dip in the charge susceptibility. Finally, for U\ = — U2, mag- 
netic and pairing moments are formed at the same temperature 
and the variations in Xspin an d Xcharge are hardly visible. 



VIII. SPECTRAL DENSITIES AND CONDUCTANCE 




Using NRG one can calculate the frequency dependence 
of the single-particle spectral densities 6 1,62 ' 63 ' 64 , which deter- 
mine the conductance through a deformable molecule with 
spin and charge degrees of freedom. We make a simplifying 
assumption that the molecule described by 2IAM is symmet- 
rically embedded between two electrodes, so that the conduc- 
tion channel in the model corresponds to the symmetric com- 
bination of electrons from both electrodes 65 . The coupling is 
then proportionate and we may use the Meir-Wingreen for- 
mula for the conductance 6 ^ 



G = Gc 



|^-Im[Tr{rG r }]^, 

ULU 



(24) 



where Go = 2e 2 /h is the conductance quantum, / the Fermi- 
Dirac distribution function, T the hybridization matrix and G r 
the retarded Green's function matrix. The components of the 
hybridization matrix are 



(25) 



where V n is the hopping amplitude from the impurity orbital 
n to the conduction band. In a simplified model we assume 
a constant density of states p and an energy-independent hy- 
bridization strength V = irpa\V\ 2 which is the same for both 
orbitals. All components of the hybridization matrix are then 
the same: Ty = T. Equation d24b simplifies to 



G = G 



-^WVlmf-i G [ 



duj. (26) 



The quantity in the parenthesis is related to the spectral den- 



sity matrix for both orbitals, Aj 



-l/(27r)Im(Gr.+G 



It 

can be computed using standard NRG techniques from matrix 
elements of the creation operators using the following spectral 
decompositions: 



A i,j (w > 0) 



- ^Re[((mKt|n ))*(m|4|n ) 
x S(u - E m ), 



A iJ (w<0) = ~^Be[((m \4\n) 

mo ,n 

x 5(lu + E n ), 



{m a \d]\n) 



(27) 



Figure 8: (Color online) Diagonal (Ai, A2) and out-of-diagonal 
spectral density (A12) of the two impurity Anderson model for differ- 
ent asymmetry parameters x. We plot only the positive frequencies. 
Due to the particle-hole symmetry, we have Ai(uj) = Ai(-ui) and 
A 12 {io) = — A 2 i(— ijj). 



where Z is the spectral sum Z = Ti-(e~@ H ), indices mo, no 
with subscript run over (eventually degenerate) ground 
states and indices m, n without a subscript over all states. 
Delta functions need to be appropriately broadened 6 ^. Note 
that there is a sum rule 



Aij{uj)dLo = 5ij, 



(28) 



which follows from the fermionic anti-commutation relation 



We are particularly interested in the symmetrized and nor- 
malized spectral density function g(u) defined by 



9 {uj) = TzT^Ai^u)) 



(29) 



where i, j = 1,2. This quantity appears in the final expression 
for the conductance: 



G = G Q 



g(u>)du>. 



(30) 



In a simple approximation, the temperature dependence of the 
conductance through the quantum dots can be deduced from 
the energy dependence of the function g(uj). 

We computed the spectral densities Ai(u>) = Ah(lj) and 
the out-of-diagonal spectral density Ai2(cj) for a constant T 
and a range of U\, U2, defined by Ui = U(l + x) and U2 — 
—U (1 — x) with constant U, see Fig. [8] For x = 0, the Kondo 
temperatures are the same. 

In the generic case x ^ (assuming x > 0, so that 
Tf ( > Tft), the spectral density of the first impurity increases 
for ui < Tf ( , but then it drops to zero for uj < T^, see the case 

of 2; = 1/5. We thus obtain a peak centered at u> ~ \Jt^-T^. 
This is at first surprising, since spin Kondo effect is usually as- 
sociated with an Abrikosov-Suhl resonance at zero frequency 
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in the spectral density and, due to the spin-charge decoupling, 
one would expect no influence of the charge Kondo effect in- 
duced by the second impurity. Indeed, at low temperatures the 
thermodynamic functions related to spin are universal func- 
tions of Tj| only; for example, x S pin in no way depends on 
the charge Kondo temperature Tjf . There is clearly no such 
universality in the dynamic quantities. Spectral density is re- 
lated to the injection of a physical fermion in the system; the 
fermion carries both spin and isospin degrees of freedom and 
is therefore sensitive to both sectors. 

For relatively large x = 1/5, the spectral density on the sec- 
ond impurity, ^(w), has a small bulge at Tj| and then it drops 
to zero. We emphasize that no Abrikosov-Suhl resonance ap- 
pears in A 2 (u}) below T% as might be expected. For x tending 
towards 0, the bulge in A2 (oj) develops into a peak of the same 
form as the one found in A\ (u>). The out-of-diagonal spectral 
density Ai2(uj) also builds up a small peak at to ~ Tk as 
x — > and goes to zero at low frequencies for all x. These re- 
sults indicate that the temperature dependence of conductance 
is expected to be non-mono tonic: it begins to increase at the 
higher of the two Kondo temperatures, but then it decreases 
when the lower Kondo temperatures is reached. The system 
does not conduct at zero temperature, g(ui — > 0) =0. This 
results is consistent with the zero temperature conduction of 
Fermi liquids as determined by the quasiparticle phase shift 
S qp = 0: 

G(T = 0) = G sin 2 6 qp = 0. (31) 

The effect of the gate voltage S is shown in Fig. [9] for the 
case of widely separated spin and charge Kondo temperatures, 
T£ <C Tj| (see also Fig. [7] where the thermodynamic prop- 
erties are shown for the same parameters). For 8 > T£, the 
isospin fluctuations on the second orbital are rapidly quenched 
and the isospin degree of freedom is fully polarized [see 
A2(u})]. The low-frequency peak in the spectral density on 
the first orbital is restored at the same time. It should be noted 
that the out-of-diagonal spectral density A\2{u>) is small for 
all 6. With increasing S, the conductance thus rises from zero 
toward the unitary limit 

G(<5»T£) ~G . (32) 

We consider the significance of such behavior in the conclu- 
sion. 



IX. TWO-IMPURITY ANDERSON-HOLSTEIN MODEL 

In previous sections, we performed calculations using the 
2IAM with a negative U 2 parameter that we considered as 
a given constant. In this section we first study how an ef- 
fective negative Us emerges in the presence of the charge- 
phonon coupling (Anderson-Holstein model) and then extend 
this model to the two impurity case with one Anderson im- 
purity (spin degree of freedom) and one Holstein impurity 
(isospin degree of freedom). 

The Anderson-Holstein model is an extension of the Hol- 
stein model in that it considers electrons with spin and (op- 
tionally) a finite on-site electron-electron repulsion in the 



impurity 14 , see Fig.Q]:. The Hamiltonian is 

H = Yj £kC L Ck ° + E ^ ( d Uka + H.c.) 

ka ka (33) 

+ U{n - l) 2 + g{a) + a)(n - 1) + u a^a 

where d is the electron annihilation operator on the impu- 
rity and a is the annihilation operator for a phonon of fre- 
quency ljq which couples to the charge n = Yla^t^ w i*h 
electron-phonon (e-ph) coupling strength g. The phonon does 
not break the particle-hole symmetry, however the full isospin 
SU(2)i sosp i n symmetry is reduced to U(l) c har g e conservation 
of charge symmetry. 

In general, the phonon renormalizes both the e-e interaction 
U, which becomes frequency dependent, and the hybridiza- 
tion strength V. In the anti-adiabatic limit loq ^S> U, the ef- 
fective e-e interaction again becomes instantaneous and the 
Anderson-Holstein model maps to the usual single impurity 
Anderson model (SIAM) witbii 

U cS = U - 2.g 2 /w . (34) 

For large enough g, U e s becomes negative and the charge 
Kondo effect is expected. In this limit (ojq ^> U), the effec- 
tive exchange interaction for the isospin degrees of freedom is 
isotropic 22 and the SU(2)i sosp j n symmetry is largely restored 
at energies much lower than gjq. 

The reversal of roles of spin and charge fluctuations in the 
Anderson-Holstein model was first deduced from dynamic 
spin and charge susceptibilities at zero temperature in Re f. [Til 
In Fig. [TO] we demonstrate this reversal in the temperature 
dependent thermodynamic spin and charge susceptibilities 
where the disappearance of the magnetic moment and emer- 
gence of the pairing moment is explicit. 

When the Coulomb repulsion in an orbital is small, U — > 0, 
but the coupling to a phonon mode is substantial, we have 
[/ ff = —2g 2 /ojQ. This holds exactly in this case (i.e. the 
condition ojq ^> U is trivially met). If, however, U and ljq in 
the same orbital are of the same magnitude, or U ^> u>o, the 
appropriate mapping is to the anisotropic Kondo mode l 15 ' 28 . 
In such systems the charge Kondo temperature is generally 
strongly attenuated, therefore in this work we do not discuss 
the case U > u>q further. 

As an illustration, we now consider the case of two active 
orbitals, the first with U\ > and the second with U2 = 
and a phonon mode, as schematically depicted in Fig.QJi. The 
Hamiltonian of this two-impurity Anderson-Holstein model 
in the particle-hole symmetric point is 

H = E £ * C L C *» 
ka 

+ f/i(m - l) 2 + .g(a t + a)(n 2 - l)+u a i a 

+ E Vk { d l ck ° + H - c ) + E Vk (4,<*a + H - c -) ■ 

ka ka 

(35) 

The presence of the phonon degrees of freedom greatly in- 
creases the degeneracy and thus the number of states that need 
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Figure 10: (Color online) Impurity susceptibilities and entropy for 
the single-impurity Anderson-Holstein model for increasing e-ph in- 
teraction strength g. U/D = 0.01, F/U = 0.04, lo /D = 0.1. 
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Figure 1 1 : (Color online) Impurity susceptibilities and entropy for 
the two-impurity model with one Holstein impurity with ujq/D = 
0.1 and varying e-ph interaction strength g, and one Anderson impu- 
rity with U/D = 0.01. Both impurities are coupled with T/D = 
0.0004 to the conduction band. 
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to be kept in the truncation step of the NRG . We used 
discretization parameter A = 3, the discretization scheme 
of Ref. 52, and we took into account the SU(2) sp i n and 
U(l) c hargo symmetries. We used f3 = 0.75 47 ' 48 . We kept up to 
5000 states (not taking into account the degeneracies) or the 
states with energy below 15T/v, whichever number is lower. 
In addition, since the eigenvalues in NRG are clustered, we 
took care not to truncate in the middle of a cluster, so that we 
do not introduce systematic errors. This approach gives con- 
verged results which agree for several values of the parameter 
z used. 

For finite g, a pairing moment is induced in the first orbital 
and a magnetic moment in the second: we then expect a simul- 
taneous spin and charge Kondo effect, as in the two-impurity 
Anderson model. This is confirmed by the results for the ther- 
modynamic properties in Fig. [TT] For small g/ujQ = 0.1, 
U e ff is small and hence the effective isospin exchange con- 
stant Jc is large: the magnetic exchange constant J$ is then 
renormalized and the spin Kondo temperature is reduced, as 
described in Sec. [HI] As g is increased, the charge Kondo tem- 
perature rapidly decreases, while the spin Kondo temperature 
returns to its unrenormalized value. At a particular value of 
g/uj Q = = 0.2236, the spin and charge Kondo 

temperatures are the same. 

X. CONCLUSION AND DISCUSSION 

We have studied the degree of the spin-charge separation in 
single-channel quantum impurity models where the impurity 
carries both the spin and isospin degrees of freedom. Using 
numerical renormalization group calculations we have con- 
firmed the well known fact that in the continuum limit, i.e. 
when the band-width D is much larger than all other relevant 
scales in the problem, spin and isospin degrees of freedom of a 
single conduction band behave as two decoupled and indepen- 
dent spin-1/2 SU(2) degrees of freedom subject only to glu- 
ing conditions. This implies that the spin and charge Kondo 
effects can coexist, the demonstration of which is the main 
result of our work. When D is comparable to other scales, 
lattice effects become important since a single orbital cannot 
sustain a spin and isospin moment at the same time. The sep- 
arate screening of spin and isospin degrees of freedom exists 
even in this limit, but the corresponding Kondo temperatures 



are reduced compared to the single-impurity results due to the 
renormalization of the Kondo exchange constants. 

We have explored the model of a deformable molecule 
with two active orbitals, one carrying a local magnetic mo- 
ment (described by the Anderson model) and one coupled to 
a strong phonon mode which induces local pairing moment 
(described by the Holstein model). We have shown that in the 
generic case the conductance through such a molecule at the 
particle-hole symmetric case has non-uniform temperature de- 
pendence: it rises at the higher of the two Kondo temperatures, 
but then drops to zero below the lower Kondo temperature. 
By changing the gate voltage from its value at the particle- 
hole symmetric point, the charge Kondo effect is quenched 
and the conductance increases back to the unitary limit if the 
spin Kondo temperature is much higher than the charge Kondo 
temperature, Tj| s> Tjf . Conversely, the magnetic field can 
be used to quench the spin Kondo effect and the conductance 
attains the unitary limit if ^> Tj|. Detection of such be- 
havior could serve as an experimental probe of the simultane- 
ous spin and charge Kondo effect. The spin-charge Kondo ef- 
fect occurs only around the particle-hole symmetric point (i.e. 
for precisely tuned gate voltage) and zero magnetic field. In 
the case of widely separated charge and spin Kondo tempera- 
tures, it would thus appear as a small region of zero or reduced 
conductance within a wider region of the high conductance 
"Kondo plateau" when plotted as a function of the gate volt- 
age and the magnetic field. In the case of similar charge and 
spin Kondo temperatures, it would be difficult to distinguish 
(at constant low temperature) the zero conductance due to si- 
multaneous charge and spin Kondo effects from the Coulomb- 
blockade valley. The temperature dependence of the conduc- 
tance is, however, different. In the case of Coulomb blockade, 
conductance monotonously decreases below the temperature 
scale of charge excitations. In the case of simultaneous charge 
and spin Kondo effects, the decrease due to Coulomb block- 
ade is followed by increasing conductance that peaks at the 
Kondo temperature before it goes to zero. 
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